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Spin splitting in the integer quantum Hall effect is investigated for a series of AbjGai-^As/GaAs 
heterojunctions and quantum wells. Magnetoresistance measurements are performed at mK tem- 
perature to characterize the electronic density of states and estimate the strength of many body 
interactions. A simple model with no free parameters correctly predicts the magnetic field required 
to observe spin splitting confirming that the appearance of spin splitting is a result of a competition 
between the disorder induced energy cost of flipping spins and the exchange energy gain associated 
with the polarized state. In this model, the single particle Zeeman energy plays no role, so that 
the appearance of this quantum Hall ferromagnet in the highest occupied Landau level can also be 
thought of as a magnetic field induced Stoner transition. 

PACS numbers: 73.43.Qt, 73.43.Nq 



I. INTRODUCTION 

The integer quantum Hall effect (QHE), 1 - observed in 
the magnetotransport properties of a two-dimensional 
electron gas (2-DEG) can be understood within the 
framework of a single electron picture. The magnetic 
field quantizes the orbital motion of electrons, so that the 
density of states consists of discrete, well defined Lan- 
dau levels, giving rise to the even integer QHE. When 
the Zeeman energy is included, the spin degeneracy of 
each Landau level is lifted, giving rise to the odd inte- 
ger QHE. While the single particle picture works well for 
the even integer QHE, this is not the case for the odd 
integer QHE, where it misses most, if not all, of the es- 
sential physics. Following the pioneering theoretical work 
of Ando and Uemura^ the important contribution of ex- 
change interactions to the spin gap is experimentally well 
established^^ 

Spectacular manifestations of many body effects have 
been extensively reported for low filling factors, notably 
at filling factor v = 1, where the ground state is an itiner- 
ant quantum Hall ferromagnet, 6 with an excitation spec- 



trum which is either a spin wavai or a spin texture exci- 
tation (Skyrmion).- The appearance of a quantum Hall 
ferromagnet, is analogous to the Stoner transition used 
to describe ferromagnetism in metals^ In contrast, at 
high filling factors, in the limit of low magnetic field, the 
Shubnikov de Haas oscillations reveal that the spin up 
and down levels are degenerate, indicating that exchange 
interactions are not sufficiently strong to induce a ferro- 
magnetic order within the highest occupied Landau level. 
In this work we investigate the transition between these 
two limiting behaviors. This transition is commonly re- 
ferred to as the appearance of spin-splitting. We show 
that this phenomenon can be understood by adapting 
the Stoner condition to take into account the modified 
2D density of states in a magnetic field. 

In the most extensively investigated system, namely 
Al^Gai^As/GaAs heterostructures, the single particle 
Zeeman energy is rather weak (~ 0.3 K/T) owing to the 
small effective electronic g-factor of GaAs (g* = —0.44 
whereas the experimental spin gaps evaluated, using 
for example thermal activation, are the order of ~ 6 K/T 
for the case of fully developed spin polarization at odd 
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filling factor. This demonstrates the importance of the 
exchange enhancement of the spin gap, which is as a re- 
sult often written as g* x HBB, where \xb is the Bohr mag- 
neton and g* x is an exchange-enhanced effective g-factor 
first proposed in Ref.^lj]. The exchange-enhanced spin 
gap depends on the spin polarization of the 2-DEG, 23 
and therefore, this description is actually restricted to 
spin polarized odd filling factors (see for example the 
theory of Aleiner and Glazman) When decreasing the 
magnetic field, experiments suggests that the enhanced 
spin gap collapses, with odd integer filling factor minima 
suddenly disappearing from the longitudinal resistivity 
(p X x(B)) at low fields. 

A theoretical explanation for the collapse of spin split- 
ting has been proposed by Fogler and Shklovskiiji 3 . who 
predicted a second order phase transition in which ex- 
change interactions are destroyed by disorder in the same 
way as temperature destroys ferromagnetism in the mean 
field theory. More precisely, when the spin gap attains 
the same order of magnitude as the disorder-induced Lan- 
dau level broadening, the spin polarization at odd filling 
factor is reduced. This leads to the destruction of the 
polarization (exchange) part of the enhanced spin gap 
which can be qualitatively understood within the early 
description of Ando and Uemura: 2 Experimentally, the 
evolution of Fogler and Shklovskii's order parameter du, 
which is the filling factor difference between two consec- 
utive p xx maxima related to spin up and down Landau 
levels, shows a universal behavior versus both tempera- 
ture and electron density after judicious re-scalingi£ii2ii£ 
Our experiments also confirms a collapse of the exchange- 
enhanced spin gap with decreasing magnetic field (Sec- 
tion WM- 

For a quantitative experimental description, param- 
eters are needed to accurately describe the competition 
between exchange interactions and disorder. This is what 
is done in Ref. |5| in which disorder is extracted from an 
analysis of the low field Shubnikov de Haas oscillations, 
and exchange is estimated from transport measurements 
of the spin gap at higher fields, when spin splitting is fully 
resolved. From this an empirical critical filling factor is 
extracted for the onset of the collapse of spin splitting. 

Here we propose to tackle the problem from an intu- 
itively different point of view. Instead of starting from a 
spin polarized situation in which the spin gap is open as 
when discussing the collapse of spin splitting, we propose 
to tackle the phase transition from the other side, by con- 
sidering what happens to an unpolarized (paramagnetic) 
2-DEG, corresponding to the Su = phase in the Fogler 
and Shklovskii approach, when increasing magnetic field. 

Spin-splitting can then be seen as a transition from a 
paramagnetic state to a ferromagnetic state within a sin- 
gle Landau level at odd filling factors. For this transition 
to occur we have to compare the energy cost, due to the 
disorder induced Landau level broadening, of populating 
higher energy levels by flipping spins (inversely propor- 
tional to the density of states at Fermi level, D(Ep)) 
and the exchange energy gain associated with the po- 



larized state. At zero magnetic field this condition is 
nothing other than the well-known Stoner condition for 
ferromagnetism in metals^ In a magnetic field, D(Ep) 
increases as B/T, where T characterizes the Landau level 
broadening. Increasing magnetic field increases the den- 
sity of states (Landau level degeneracy), thus lowering 
the energy cost for flipping spins sufficiently to induce a 
ferromagnetic state. For a quantitative comparison with 
experiment, estimations of the Landau level broadening 
r and the strength of the exchange energy are required. 

The rest of this paper is organized as follows. In Sec- 
tion [H] a simple model for the appearance of spin split- 
ting, based on the Stoner condition for ferromagnetism 
is presented. The Landau level broadening is evaluated 
in Section IlIII from an analysis of the low magnetic field 
Shubnikov de Haas oscillations. The exchange energy is 
derived in Section Hvl using the theoretical calculation of 
the electronic spin susceptibility of Attaccalite et alJ&. 
The calculated spin susceptibility and exchange parame- 
ter are shown to be in good agreement with experiment. 
Finally, in Section [3 the Landau level broadening and 
exchange energy are used to predict the critical magnetic 
field for the appearance of spin splitting. A comparison 
with the experimental data, shows good agreement for a 
number of samples covering a wide range of densities and 
mobilities. 



II. MODEL FOR THE APPEARANCE OF SPIN 
SPLITTING 

We propose here that exchange interactions drive a 
transition from a paramagnetic state, corresponding to 
non spin-split Shubnikov de Haas oscillations, to a fer- 
romagnetic state, corresponding to spin resolved Shub- 
nikov de Haas. Strictly speaking, the transition is driven 
by the density of states, via the Landau level degeneracy 
eB/h, rather than the exchange energy, which is actually 
magnetic field independent, to a good approximation, at 
high filling factors (see Section ITV _B|) . The system is de- 
scribed considering the highest occupied Landau level at 
odd filling factor, so that the starting situation is when 
the system is unpolarized and the Fermi level lies in the 
center of the degenerate spin up and down sub levels of a 
given Landau level. This is an approximation, since there 
can be a small residual spin polarization due to the Zee- 
man splitting. However, in GaAs, at low magnetic fields, 
the Zeeman energy is small compared to both Landau 
level broadening and exchange energy, so it is a good ap- 
proximation to consider the limit of zero Zeeman energy. 
Within this framework the appearance of a ferromagnetic 
state requires that the energy cost of populating higher 
energy levels by flipping spin, should be less than the gain 
in exchange energy, which therefore stabilizes the newly 
polarized state. This situation is depicted schematically 
in FigUl 

If we consider flipping spins to fill the spin up Landau 
level to an energy u above the initial Fermi energy, the 
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FIG. 1: Schematic representation of the process for spin split- 
ting. 



energy cost can be expressed as, 

E fUp = 2 f D(E)E dE. (1) 
Jo 

To simplify, we have defined the zero of energy such that 
Ep = 0. It is possible to obtain a simple expression 
if we assume that the density of states, D(E), in this 
expression can be approximated by its value at the Fermi 
level. Proceeding in this way, we implicitly suppose that 
the Landau level is rectangular shaped, which should be 
relevant only for small values of u. If the temperature 
is zero and to = 2 / " D(E)dE w 2D(E F ) u is the spin 
polarization created by this process within Landau level 
N, then Eq.JTJ can be written as, 



E 
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4D(E F ) 



(2) 



The factor 4 in Eq.J2J simply arises from the definition 
of to. This expression confirms the intuitive idea that 
the larger the density of state is, the smaller the energy 
cost of flipping spins will be, simply because more states 
are available at lower energies. Calculations assuming a 
Gaussian form for the Landau level show this approxima- 
tion differs from the exact result of Eq.QJ by less than 5 
%, provided to < 0.425. 

At the same time this spin flip process leads to an 
increased population of the spin up Landau level, giving 
rise to an exchange energy gain which can be written, 



E x = —m 



(3) 



where Xn is the exchange energy between 2 spins in Lan- 
dau level N. The quadratic dependence on to results 
from the many body nature of this interaction. Follow- 
ing this model in which we only have two competing en- 
ergies, it will be energetically favorable for spins to flip 
when, 
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When the Fermi level lies the center of a Landau level, 
D(E F ) increases as B/T, due to the eB/h Landau level 
degeneracy and the disorder-induced broadening T. The 
underlying idea is that when increasing the magnetic field 
we reach a sufficient density of states at the Fermi level to 
reduce the spin flip energy cost, so that the spin system 
eventually becomes "ferromagnetic" . 

To check the validity of such a model, one must eval- 
uate quantitatively how the density of states varies with 
magnetic field. For this a description of the shape and 
width of the Landau levels is required. 



III. DENSITY OF STATES IN A MAGNETIC 
FIELD 

A. Low field Shubnikov de Haas analysis 

1. Formalism 

In a magnetic field, the diagonal conductivity a xx of 
the 2-DEG can be written, 



5 6 Ttr 



771* (1 



(5) 



where n 8 is the electron density, T tr the transport scat- 
tering time, to* the electron effective mass, and lo c the 
cyclotron frequency. 

In the limit of a zero temperature and if the magnetic 
field is large enough to satisfy Lu c T tr 3> 1 , this simplifies 
to give, 



e 2 D{E F 

m*T tr UJr 



(6) 



For high mobility samples, the Hall resistivity is usually 
much larger than the longitudinal resistivity provided 
u c Ttr 3> 1, so that the 2-dimensional resistivity tensor 
gives p xx « a xx p xy . In the non spin-split low magnetic 
field region, p xy is approximatively linear in B, and from 
Eq.@ the field dependance of p xx can be written, 



Pxx oc —D(E F ). 

Ttr 



(7) 



At low magnetic field, the variation Ap of the longitudi- 
nal resistivity around the zero field resistivity (po) can be 
calculated by expanding the quantized density of states 
as a Fourier series. The amplitude of the fundamental 
oscillatory term gives, in the zero temperature limit, the 
following expression for Ap/ p^, 



Ap 
Po 



Ae z 



(8) 



This expression, which is a good approximation pro- 
vided that Ap/po < 1, is obtained for the particular 
case of Lorentzian Landau level and is nothing more 
than the disorder damping term in the Lifshitz Kosevitch 
formula. 1 '' Here T q is the quantum life time related to the 
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width of the Lorentzian Landau level by T q = h/2T l ■ A 
is a constant whose value depends on the power relation 
between the resistivity and the density of state (see be- 
low). 

The same calculation for a Gaussian Landau level 
(defined here with a full width at half maximum of 
2^1n(2)r G ) gives, 



Ap 

Po 



Ae 



(9) 



The expression is similar to that for a Lorentzian Lan- 
dau level except that the B dependence in the exponen- 
tial function is in 1/ B 2 if the Landau level width is field 
independent. 

Plotting the experimental Ap/ po on a logarithmic scale 
versus 1/B or l/B 2 can give information concerning the 
shape of the Landau level. A linear behavior of ln(Ap/po) 
when plotted versus 1/B (so-called Dingle plot) indicates 
Lorentzian Landau levels^, whereas a linear behavior 
when plotted versus 1/B 2 suggests Gaussian Landau lev- 
els. In both cases the slope of these curves gives the 
(field-independent) width of these levels. This theoret- 
ical treatment has however to be confirmed by check- 
ing the value of the intercept at 1/B = 0, noted here 
as A. An intercept of A = 2 is theoretically expected 
when the resistivity is directly proportional to D(Ep). 
If however, the p xx dependence on D(Ep) is quadratic 
(i.e. l/r fr oc D(Ep) in Eq.Q), the Fourier expansion 
exhibits an additional factor 2, which leads to an inter- 
cept of A = 4. These issues are discussed in more detail 
in the seminal paper by Coleridge^ 



2. Low temperature magnetotransport 

We have performed mK temperature magnetotrans- 
port measurements on different hetero junctions(HJ) and 
quantum wells (QW), using a standard low frequency 
lock-in technique. Representative longitudinal resistivity 
(Pxx) versus magnetic field data obtained for one of the 
samples is shown in Fig|3 The variation of the longitu- 
dinal resistivity around the zero field value (Ap/po) has 
been directly extracted from the low field Shubnikov de 
Haas data at 50mK. (Ap/p ) is plotted on a logarith- 
mic scale versus (1/B 2 ) in Fig[3fa), and, for comparison, 
versus (1/B) in Fig|2Ib), for all the samples investigated. 

At T = 50 mK, the smearing imposed by the Fermi- 
Dirac statistics is small enough to ensure that Ap can 
be analyzed using a zero temperature formalism. In the 
samples studied, the Dingle plots are found to decrease 
as 1/B 2 rather than as 1/B (see Fig[3J). While the linear 
fits are reasonable in both cases, they are slightly better 
in the 1/B 2 plots (FigEJa)). The data in the 1/B plots 
has a noticeable curvature (Fig^b)). In addition, the 
intercepts in the 1/B 2 plots (1 < A < 2.7) are close to 
one of the theoretically expected values, A = 2. This is 
not the case for the 1/B plots, with intercepts 7 < A < 




em 



FIG. 2: Magneto resistance trace for sample NU 1783a. 8v 
(full circles) as a function of magnetic field. Dotted line is 
the extrapolation of 8v used to determine B ss as explained 
in Section [V] Vertical solid line delimits the magnetic field 
region for which Ap/ po < 1. Dash-dotted straight line em- 
phasizes the linear behavior of p^ k (B) at higher magnetic 
fields. 



21, far from the theoretically expected value of either 
A = 2 or A = 4. 

This suggests, that for our samples, the density of 
states is better described using a Gaussian form with 
field-independent width for the Landau levels. Within 
this model, the intercept A ~ 2, implies that p xx is di- 
rectly proportional to D(Ep). We have also checked that 
an identical result is obtained when filtering the raw data 
with a band pass centered at the fundamental Shubnikov 
de Haas frequency. The 1/B 2 decrease (a Gaussian form 
for the Landau levels) therefore provides a better descrip- 
tion of our data. From Eq.Q this result implies, as pro- 
posed in Ref. [2(j that the transport scattering time does 
not depend on D(Ep). 

A complementary study to the low field Shubnikov 
de Haas analysis can be performed by analyzing the 
values of the resistivity peaks, (P xx ak (B), in the regime 
Ap/ po > 1,22*21 In this regime, we observe a linear in- 
crease of p v xx k (B) (visible in the p xx trace for sample 
NU1783a in Fig©. When Ap/p > 1, Eq.© and Eq.@ 
are no longer valid, because the deviation of the density 
of states from the zero field value becomes too large to 
be described by the fundamental oscillatory term in the 
Fourier series. However, Eq.(|7|) is still valid and thus the 
field dependence of the resistivity peaks can provide in- 
formation on -^-D(Ep). The less stringent condition to 

account for the observed linear dependence of p x \ e x ak (B) 1 
is as proposed in Ref.[2fJ, that the inverse of the trans- 
port relaxation time and the Landau level broadening 
have the same mag netic field dependence. For the sample 
studied in Ref. 20], this dependance was in square root 
of B, while here the low field Shubnikov de Haas (Dingle) 
analysis suggests both parameters are field-independent. 
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FIG. 3: Dingle plots, (a) Ap/po versus 1/B 2 and (b) Ap/p 
versus 1/B. Note log scale for vertical axes. Straight lines 
are linear fits to the data points in the region Ap/ po < I 
delimited by the dotted horizontal dotted line. 



The latter case is fully compatible with the observed lin- 
ear dependence of p x l e x ak {B). 

The theory of Raik and Shabazyar. 22 on the Lan- 
dau level shape in the presence of long range scattering, 
suggests Gaussian Landau levels with field independent 
width in the regime where a corr 3> R cy d (where a corr 
is the correlation length of the disorder potential, and 
Rcyd the cyclotron radius) and with a square root of B 
dependence when R cyc i 3> a corr . Our results are con- 
sistent with a field-independent broadening of the Gaus- 
sian Landau levels suggesting that a corr ^> R cyc i in our 
samples. While this is an interesting point discussed in 
Refs. [l9ll2^ |. it is not the direct focus of the present work. 



B. Estimation of the Landau level width 

To evaluate the D(Ep) appearing in Eq.QJ, one needs 
an estimation of the total Landau level width T, which in- 
clude all states, both delocalized and localized. Estimat- 
ing r from the low field Shubnikov de Haas oscillations 
is a priori not possible because in this region delocalized 
states from different Landau levels overlap, preventing 
the resistance from going to zero. However, at slightly 
higher fields, in the region just before spin-splitting ap- 
pears, the proximity to the first zero resistance states at 
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FIG. 4: Schematic representation of the density of states in a 
magnetic field used in a simple model to estimate V in Section 

ITTTI 



even filling factors can be used to estimate T. The width 
of the p X x minima (zero resistance states) depends on 
the ratio Tdi/T, where Tdi characterizes the width of the 
extended (delocalized) states. 

Starting from the formalism presented in section llll Al 
we construct a simple model in the zero temperature 
limit, assuming that the density of states can be mod- 
eled as a set of Gaussian functions of width T, with sub- 
Gaussians delimiting delocalized states near the center of 
the Landau level, of width Tdi ■ This is schematically de- 
picted in Fig0J Tdi = ^g, is determined from the slope of 
the ln(Ap//?o) versus 1/B 2 plots as described in section 
IIII Al The values of for all the samples investigated 
are summarized in table [I] p xx can then be calculated 
using p xx oc D(Ep). Knowing Tdi, it is possible to de- 
termine r, for each sample, by fitting to the p xx {B) data 
in the magnetic field region just before spin-splitting ap- 
pears. Proceeding in this way we have estimated T for the 
different samples studied (see table QJ. Experimentally a 
Gaussian form for the density of delocalized states gives 
by far the best fit to our p xx (B) data. We have also 
tried to fit assuming the region of delocalized states is 
delimited with a sharp cut-off (mobility edge), as often 
proposed in the literature. While this model is unable to 
correctly reproduce the shape of p xx (B) we can never- 
theless roughly estimate T. The results are comparable 
to those obtained with our Gaussian model, suggesting 
that the Landau level width T extracted here is to some 
extent independent of the exact form used for the density 
of delocalized states. 



IV. ESTIMATION OF THE EXCHANGE 
INTERACTION 

Analyzing the spin split Shubnikov de Haas oscillations 
we can extract parameters which characterize exchange, 
such as, for example, the enhanced-effective g-factors^ 
However, in this case the spin gap already exists, and, 
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Sample 


Structure 


n s 


A* 


r M (K) 


r (K) 


B ss (T) 


9ex 






(10 n cm- 2 ) 


(10 6 cm 2 /Vs) 










F1201 


QW 


8.83 


1.0 


1.23 ±0.05 


2.3 ±0.1 


0.93 ±0.07 


5.3 ± 1.5 


F1200 


QW 


7.55 


1.0 


1.35 ±0.05 


1.9 ±0.1 


0.72 ±0.10 


5.7 ±0.4 


NRC1707 


HJ 


1.64 


5.8 


1.10 ±0.05 


1.5 ±0.1 


0.26 ±0.05 


13.1 ± 1.0 


NU 1783a 


HJ 


2.10 


1.5 


1.48 ±0.05 


2.2 ±0.1 


0.43 ± 0.06 


13.0 ± 1.0 


NU 535 


QW 


2.13 


0.11 


2.00 ±0.05 


3.1 ±0.2 


0.74 ±0.12 


9.5 ± 1.0 


NU 1783b 


HJ 


1.77 


1.8 


1.53 ±0.05 


2.3 ±0.2 


0.32 ±0.11 


13.7 ±0.7 


N 178 


HJ 


1.32 


0.32 


2.50 ± 0.05 




0.54 ±0.2 


16.0 ± 1.5 



TABLE I: Samples parameters: structure, electron density (n s ), mobility at 50mK (/x), width of delocalized states (Tdi), total 
Landau level width (r), magnetic field for the appearance of spin-splitting (B ss ) discussed in Section[V]and exchange-enhanced 
g- factor (gt x ) discussed in Section fl V Bl For N 178 a reliable estimation of F was not possible. 



the system is therefore already polarized. To describe an 
unpolarized 2-DEG in a low field situation, it is better 
to estimate exchange via the spin susceptibility of the 
electronic system, defined as the response to an external 
magnetic field in terms of spin magnetization. This is 
a good probe of exchange interactions between electrons 
since exchange favors spin alignment and thus enhances 
the spin susceptibility. We will see in section IV. B. that 
these two different approaches to estimate the exchange 
(enhanced-effective g-factors and spin susceptibility) are 
fully consistent. 



A. spin susceptibility 

1. Theoretical approach 

At zero magnetic field, the nature of the 2-DEG ground 
state is strongly dependent on the electron density, and 
thus often described as a function of the so-called density 
parameter r s , corresponding to the ratio of the Coulomb 
energy to the kinetic energy. In 2D, r s scales as l/^/riJ. 
At high r s , in the low density regime, the ground state 
remains a topic of controversy, with a possible manifes- 
tation of a metal to insulator transition (MIT) observed 
around r s ~ 10^ At low r s however, and especially for 
the density range investigated here (r s < 2), the density 
is sufficiently high to ensure that the zero field ground 
state is a paramagnetic liquid. In this case, the para- 
magnetic relative spin susceptibility x/xo of the 2-DEG 
can be derived from calculations of the B = ground 
state energy as a function of density and polarization, 
for example using quantum Monte Carlo simulations (see 
e.g. Tanatar and Ceperley )j24 Here, x/xo corresponds 
to the ratio of the "real" spin susceptibility of interact- 
ing electrons, x, to the non-interacting Pauli value xo- 
Similar calculations have been performed more recently 
by Attaccalite et alw&. Different predictions arise from 
the different ways of estimating the total energy of the 
ground state, and in particular the correlation term. At 
low density, x/xo turns out to be a very delicate quantity 
to estimate, and the various predictions are significantly 
different. 

Fortunately, for our rather high density systems (r s < 



2), the behavior of the predicted susceptibilities is more 
robust. Indeed, predictions for small values of r s do not 
show large discrepancies because in the limit r s — > 
they all reach the "basic" Hartree-Fock approximation 
for spin susceptibility, determined using the well- known 
Hartree-Fock 2D ground state energy (E g — l/r 2 s — 
1.2004/r s , see e.g. Ref.|H). This Hartree-Fock un- 
screened susceptibility increases with increasing r g , re- 
flecting that the role of exchange interactions becomes 
more important at low density. This trend is common to 
a large majorities of predictions over a very large r s do- 
main. However, the Hartree-Fock susceptibility diverges 
for r s ~ 2.2 because it does not include the effects of 
the correlation energy and/or screening, which prevents 
x/xo from diverging. Taking into account these correc- 
tions, as done for instance in Ref.0], the r s increase is 
still present but much smoother. These different theoret- 
ical predictions are plotted for comparison in Fig|SJ 



2. Experimental approach 

There are several experimental techniques to measure 
the electron spin susceptibility of a 2-DEG (for a review 
see Ref . 26] ) . For the particular case of a 2-DEG in GaAs, 
significant results have recently been obtained with the 
tilted field method by Zhu et alSL. The principle of this 
method is to induce Landau level coincidences by tilting 
the sample in a magnetic field. When tilting the sample, 
the magnetic field perpendicular to the 2-DEG is reduced 
(see inset of Fig. EJ) . We recall that the Zeeman gap de- 
pends on the total magnetic field, whereas the cyclotron 
gap, which is an orbital effect, depends only on the field 
component perpendicular to the 2-DEG. 

Therefore, by tilting the sample, for a given filling fac- 
tor, we can increase the relative weight of the Zeeman 
gap until the point where the spin gap is equal to the cy- 
clotron gap. This is what we refer to here as the coinci- 
dence condition (note in the literature coincidence condi- 
tion often refers to the case where the spin gap is half the 
cyclotron energy). Knowing the angle 9 (defined in the 
see inset of Fig. [5| which corresponds to this situation, 
the product m* gl ilt satisfying m* g* t ut = *2m e cos(#), 
where i = 1 can be extracted. If normalized using the 
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band structure parameters (m* = 0.068m e , m e being the 
electron rest mass, and \g*\ — 0.44 ), this m*g^ iH is di- 
rectly equal to the relative spin susceptibility x/xo a t 
i = 1. The parameter g^ ilt describe the total spin gap 
assuming it can be written g^ ilt fiBB where B is the to- 
tal magnetic field consistent with the definition of the 
spin susceptibility. The enhanced g-factor can also be 
defined using the perpendicular component of the mag- 
netic field, related to the 2-dimensional (orbital) nature 
of exchanged 

Proceeding with a tunable density gated sample, Zhu 
et alml give an explicit density/polarization dependence 
of the relative susceptibility for a range of density go- 



ing from 2 x 10 9 to 4 x 10 1 



The polarization de- 



pendence is obtained from the observed i-dependence of 
x/xoj because different values of i correspond to a dif- 
ferent polarization of the system. The limit of zero po- 
larization can provide an estimation of the paramagnetic 
relative susceptibility of the 2-DEG. Their results show 
the experimental x/xo is much smaller than the Hartree- 
Fock prediction and also that the r s trend is in good 
agreement with the Attaccalite et al. calculation, al- 
though the experimental values are slightly smaller than 
predicted. Other experimental data on different systems, 
mainly 2-DEG in Si, can complete this comparison if they 
are rescaled in terms of r s (see e.g. Rcf. 26]), and show 
that this trend is also respected in the lower r s region. 

To check if Attaccalite et al. prediction is experimen- 
tally valid for GaAs in the region of interest here (r s < 2), 
we have performed tilted-field measurements on our sam- 
ples which cover a higher density range than the one 
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FIG. 5: Theoretical spin susceptibility as a function of the 
density parameter r s . Hartree- Fock approximation (dotted 
curve), Attaccalite et al. prediction (solid curve). Open 
squares are experimental values of x/xo at i = 1 determined 
by tilted-field experiments on different samples. For the high 
density sample (r s ~ 0.7) only an upper limit could be ex- 
tracted. Full circles are extrapolated (B — > 0) values of the 
paramagnetic x/Xo- 
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FIG. 6: Magnetoresistance traces for sample NU 1783b at 
different tilting angles. Angles are determined from the slope 
of the low-field Hall resistance. 



in Pi ef . [271] . Going towards high density should reduce 
spin susceptibility and thus requires even larger tilt an- 
gles, and thus higher magnetic fields. Transport mea- 
surements were performed using standard low frequency 
lock-in technique under magnetic fields up to 23 T, in 
a dilution fridge (base T w 50 mK) equipped with an 
u in-situ" rotating sample holder. 

In FigEl we plot p xx (B), measured at T = 50 mK, 
and at various tilt angles, for sample NU1783b. Note 
different traces are "naturally" shifted upwards, due to 
the effect of the increasing in plane magnetic field which 
leads to a positive magnetoresistance^ The angles were 
accurately determined from the slope of the low mag- 
netic field Hall resistance. As expected, we find that 
coincidence is reached for even larger tilt angles than in 
the Zhu et al. experiment^ This is consistent with pre- 
vious measurements in this higher density ranged For 
example in sample NU1783b coincidence is reached at 
approximately 88.7° (see FigEJ. The fact that we can 
use tilted field to tune the Zeeman energy is not incon- 
sistent with our simple model in which the single particle 
Zeeman energy is assumed to be negligible in the 9 = 0° 
configuration. The large tilt angle required, implies that 
an in plane magnetic field approximately two order of 
magnitude larger than the perpendicular magnetic field 
component is necessary to achieve Landau level coinci- 
dence. 

If we consider the situation at even filling factor, for 
example, at v = 16 indicated by the arrow in FigEI we 
see that when increasing 6 the gap at even filling factor 
initially closes, until the spin up and down sub-levels of 
adjacent Landau levels coincide, here at about 88.7 de- 
grees. Tilting further, we are able to make the levels 
cross, confirming we do indeed obtain a coincidence of 
Landau levels. We can see that the gap re-opens at even 
filling factors, and then closes at odd filling factors prob- 
ably reaching a second coincidence in which the spin gap 
is twice the cyclotron gap. The coincidence condition de- 
pends only weakly on the perpendicular magnetic field. 
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In other words, a single tilt angle 9, provides coincidence 
nearly simultaneously at all observed even filling factors. 
This means that the polarization dependance of x/xo is 
weak, more precisely sufficiently weak, that its effects are 
suppressed by Landau level broadening. 

Determining the coincidence condition from the p xx 
trace, the corresponding products m*gi at have been ex- 
tracted for two samples measured (1783b and N178). For 
the high density samples, coincidence occurs at higher 
tilt angles (9 > 88.81°), and, due to sample mobility and 
magnetic field limitations, coincidence cannot be reached. 
Nevertheless, we can extract an upper limit for the prod- 
uct m*g% ilt of sample 1200, showing that the susceptibil- 
ity is indeed slightly lower for this higher density sample. 
The values of m*g^ ut , normalized by the product of the 
band structure values m* = 0.068m e and \g*\ = 0.44, are 
plotted versus r s in Fig0(open squares). 

The paramagnetic spin susceptibility x/xo can then be 
estimated from the i=0.5, 1 (and 1.5 when it is observed) 
conditions at a given perpendicular magnetic field, as- 
suming the i-dependance of spin susceptibility is linear as 
observed in Ref . [23 ■ Such a linear extrapolation induces 
large error bars, because in the best situation we only 
have 3 values for x/Xo(i)t themselves subjected to error. 
As the tilt angle required for different coincidence condi- 
tions is roughly the same for all filling factors, extrapola- 
tion at different perpendicular magnetic field gives sim- 
ilar results. The differences are within the error bars 
induced by the linear extrapolation. While, the value of 
x/xo a t i = l is subject to less uncertainty than the ex- 
trapolated paramagnetic x/xo, rigourously it is the latter 
which is required for a comparison with theory. 

The paramagnetic x/xo ar e plotted versus r s for sam- 
ple 1783b and N178 in Fig[5] (closed circles). They are 
consistent with the Attaccalite et al. prediction, smaller 
at most by 20 %. The experimental values in Ref.|27j are 
also smaller than the theoretical prediction. However, 
the calculation of Attaccalite et al. is performed for an 
ideal 2-DEG, whereas experimentally the finite thickness 
of a real 2-DEG is expected to reduce spin susceptibility. 
Very recent finite thickness corrections to the Attaccalite 
et al. calculation have been performed by De Palo et 
alm&, showing these effects indeed lead to a reduction of 
the theoretical values, reconciling them with experiments 
(see also Ref. [30l] ) . Furthermore, spin susceptibility mea- 
surements on system similar to the one used in Ref. |27j| 
have been recently extended to the high density regime 
down to r s = 0.8^i showing good quantitative agree- 
ments with the finite thickness correction in this density 
regime 
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What is important for this work, is that the Attaccalite 
et al. calculation provide a relevant estimation of spin 
susceptibility in our samples, and therefore we can use 
this theory to derive our exchange parameter. 
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FIG. 7: Measured p xx versus magnetic field for sample 
NU1783a at T = 50 mK (solid line). Calculated p xx taking 
into account exchange interactions with a phenomenological 
spin gap g* ex psB (broken lines). 



B. Exchange parameter versus density 

We now have to draw a formal link between the rela- 
tive spin susceptibility and the exchange parameter Xn 
defined in Section |H] As x/xo is calculated for the case 
of zero magnetic field, our condition for the appearance 
of spin splitting should be applied in the limit B = 0, 
which gives X oo D(Ef)b=0 = 1- 

As this condition correspond to a zero-field "ferromag- 
netic" state, we can define, 



X_ 
Xo 



1 



1 - X 00 D{E F ) B =o 



(10) 



Here D(Ef)b=q is the zero-field density of states of the 
2-DEG, for one spin orientation, i.e. m* /2ttTi 2 . Here 
m* is taken to be equal to the accepted GaAs band edge 
value of 0.068m e . The zero field exchange parameter X^ 
is then obtained from the Attaccalite et al. spin suscep- 
tibility (x/Xo)au as a function of density, 



1 



{%)Att(n.)D(E F ) B=0 



(11) 



In order to compare this prediction with experiment, 
we have estimated the exchange-enhanced effective g- 
factors g* x . Starting from the density of states, D{E), 
as described in Section IIIII we use a simple zero tem- 
perature model, in which D{E) is described by a set of 
Gaussian Landau levels with magnetic field-independent 
width. 3 - The longitudinal resistivity p xx is calculated us- 
ing p xx cx D(Ep). We have extracted the value of g* x 
by fitting odd filling factor minima in p xx assuming that 
the Landau level broadening involved here is the same 
the one extracted from the low field Shubnikov de Haas 
oscillations. This assumption, also used in Ref.[|J, sup- 
poses that the field independence of the Landau level 
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FIG. 8: Exchange parameters as a function of the density 
parameter r s , and electron density n a (note nonlinear scale 
for the top axis ). Estimation from Ando and Uemura" (dash- 
dotted line), from Aleiner et aZi—(dashed line) and from the 
experimentally determined g* x (full circles). The dotted line 
is a linear fit to gt x {r s ). 



width, observed at low field, still holds in the magnetic 
field region where spin-splitting first appears. This field 
independent Landau level width is suggested by the 1/B 2 
dependence of ln(Ap/p ) in the Dingle plots (Fig|3J. 

As expected, a single electron picture, in which the 
spin gap is taken to be equal to the Zeeman gap, fails 
to reproduce the observed spin splitting. It is however, 
possible to reproduce odd filling factor minima using a 
magnetic field (filling factor) dependent phenomenologi- 
cal spin gap, g* x fisB, where g* x takes into account the 
exchange enhancement of the spin gap. Representative 
results can be seen FigF7| where we plot p xx {B), mea- 
sured at T = 50 mK, for one of the samples. The broken 
lines are the simulated p xx (B) calculated for each odd 
filling factor using the g* x indicated in the legend. The 
exchange-enhanced effective g-factor collapses when ap- 
proaching the non-split p xx region, consistent with the 
Fogler and Shklovskii phase transition^H reflecting the 
disorder-induced destruction of spin polarization at odd 
filling factor. 

When the spin gap is fully open, the spin polarization 
is a maximum and the value of g* x (reported in table P) 
can be compared with theory. Neglecting the single par- 
ticle Zeeman energy, we can equate the exchange gap, 
X(eB/h), to the total spin gap gl x psB. This corre- 
sponds to the fully polarized case in which a spin interacts 
with the eB/h electrons in the lower spin Landau level. 
The values of X = g% x PB{h/ e) for the different samples 
studied are plotted in FigEl (closed circles) together with 
the predicted value Xqq obtained by Eci. lTfTJl (solid line). 
As previously observed in Al x Gai^ x As/GaAs (see e.g. 
Ref.[5j), g* ex is larger in low-density samples. The esti- 
mation derived by Ea. i|ll|) from Attaccalite et al. spin 



susceptibility is in good quantitative agreement with the 
one obtained from our g* x , even though the exact func- 
tional dependence upon r s (sub-linear) is not well repro- 
duced by the data. 

For comparison, we also show in Fig|H|other theoretical 
estimations of the exchange parameter. Using the Ando 
and Uemura— exchange-enhancement of the spin gap, we 
can again write X(eB/h) — ^(e 2 kp /Aire). Here the mag- 
netic length has been replaced by 1/kp which is the rel- 
evant length scale at high filling factorsAi^. The filling 
factor v = hn s /eB and kp — ^2nn s so that in Fig|S]we 
plot X = e 2 y / 2'ir/n s /Aire (dotted-dashed line). The de- 
pendence on 1/kp rather than Ib leads to a cx r s 
density dependance. It is striking how this simple argu- 
ment describes qualitatively the density dependance of 
our enhanced g-factor g* x tx r s (as visible on the linear 
fit to g* ex (r s ) in dotted line). We have also estimated the 
exchange parameter from the a parameter describing the 
exchange gap in Ref . |l2( , which is valid for r s < 1 (dashed 
line). The rigorous calculation^ including screening, in- 
troduces a logarithmic correction which reduces the ex- 
change parameter at higher r s , improving agreement with 
experiment. 

Finally we note that the influence of disorder on the 
measured g* x seems to be weak since different samples 
with different Landau level broadening keep within the 
density trend. This is consistent with the fact that the 
effect of disorder on the spin susceptibility extracted by 
the tilted field method is known to be weakiSi To summa- 
rize, the good agreement observed between experiment 
and theory, means that the exchange parameter for our 
samples can be reliably estimated using Attaccalite et al. 
spin susceptibility. To a good approximation, the Xn 
appearing in Eq.Q can be identified with its zero field 
value, Xoo, given by Eq.lfnj. The absence of a magnetic 
field dependence can be understood, since at high filling 
factors (low B) the average electronic separation is 1/kp, 
rather than the magnetic length. 



V. CRITICAL MAGNETIC FIELD FOR SPIN 
SPLITTING 

A. Predicted field for spin splitting 

From Eq.(T4} in section [HI we can directly obtain a 
relation which explicitly gives the magnetic field B ss at 
which the spin degeneracy of a Landau level should be 
lifted, as a function of the Landau level broadening T, 
and the electron density n s , 



B, 



e Xoo(n s )' 



(12) 



where X^ is given by Ea. HllI) . We see here that B ss is 
predicted to increase with T, so that spin splitting ap- 
pears later in more disordered samples, which is also ex- 
perimentally well-established. The same T-dependance 
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is found for the condition for the collapse of spin split- 
ting when comparing the spin gap to the Landau level 
broadeningi^ii The electron density dependence of B ss 
arises via the dependence of on n s . With de- 
fined as in Eq.lJTIJl this gives for B ss a dependence close 
to .y/nj, pushing spin splitting to higher magnetic field 
for higher density samples. 



B. Comparison with experiment 

Experimentally, the spin splitting phenomenon is 
probed by the presence of a minimum appearing at odd 
filling factor in the longitudinal resistance. A method of 
extracting B ss from the experimental data is needed for 
a quantitative comparison. In Ref. [l3l ] , the critical filling 
factor (magnetic field) for collapse of spin splitting is de- 
fined as the filling factor corresponding to 5v = 0.5. Here 
5v, is the filling factor difference between two consecutive 
p xx maxima related to spin up and down Landau levels. 
This method, has been successfully applied to analyze 
the "collapse" of spin-splitting, in which the spin gap is 
compared to the disorder in order to determine when the 
spin splitting should collapsed 

In a similar manner, we have analyzed our p xx {B) data 
to extract 8v, as shown in Fig|2]for sample NU1783a. In 
our approach, in which we start from the "paramagnetic" 
state, the condition 8v = is more appropriate. As 
the condition Sv = cannot be accessed experimentally, 
we have to extrapolate from the experimental 8v{B) 1 for 
which we typically have data in the regime 0.5 < Sv < 1 
(see Fig|2J). The value of B ss determined in this manner 
is only slightly smaller than if we had used the Sv = 0.5 
criterium, and, the error implied by such a treatment is 
estimated for each sample from the width of the transi- 
tion. 

Rigourously, as our model is developed for the zero 
temperature case, what is required is the extrapolated 
value of B ss at zero temperature. However, we have per- 
formed this analysis at different temperatures (0.05 < 
T < 1.2 K) and the extracted evolution of B ss versus 
temperature shows that the difference between the ex- 
trapolated T = K value, and the 50 mK value, never 
exceeds 5%, which is well within the error bar for B ss . 

In Fig|51 B ss is plotted both versus Landau level broad- 
ening (F) and electron density (n s ), for the six differ- 
ent samples (HJ or QW), spanning a density and mo- 
bility range respectively of ~ (1.5 — 9) x 10 11 cm -2 and 
~ (0.1 — 6) x 10 6 cm 2 /Vs. The T values, which are also 
given in Table [0 have been estimated as explained in 
Section ITTll To compare experiment with our model, the 
predicted evolution of B ss calculated using Ea. ll2|) . is 
plotted as a wire mesh (3D plot) in Fig|5| The exchange 
energy X^ has been estimated using the calculated spin 
susceptibility of Attaccalite et alm& as explained in Sec- 
tion ^3 FigEl shows a good quantitative agreement be- 
tween the prediction and the experimental B ss , especially 
considering that there are no adjustable parameters in 
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FIG. 9: Magnetic field for spin splitting (B ss ) plotted versus 
Landau level broadening (r), and, density (n s ). The wire 
frame is the prediction of Eq. ll2ll . and the spheres are the 
experimentally determined B sa . The experimental error is 
indicated by the size of the spheres. 
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FIG. 10: Critical density of states at the Fermi level 
(D 3S (Ef)) for a paramagnetic-ferromagnetic transition plot- 
ted as a function of the inverse exchange energy 1/Xoo. The 
Dss(Ef) predicted as explained in the text using Eq. H2H 
(solid line) and the zero field 2D density of states (broken 
line) are also plotted for comparison. 



the model . If we focus on the B ss variation at con- 
stant disorder, the density dependence is correctly de- 
scribed. Equally, at fixed density, the disorder depen- 
dance is roughly linear in T as predicted by Ea. (|12fl . 

The good quantitative agreement gives further weight 
to the main idea of our model in which there is a compe- 
tition between the energy cost of a flipping spins, which 
increases with T, and the exchange gain of flipping spins, 
which increases with decreasing n s , because of stronger 
exchange interactions at lower densities. The assump- 
tion that spin splitting in GaAs is driven primarily by 
exchange interactions, the Zeeman splitting being only a 
correction, is also validated. 

From the experimentally determined critical magnetic 
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field, B ss , and the Landau level broadening, T, it is possi- 
ble to calculate the critical density of states at the Fermi 
level, D ss (Ef) = {eB 8B /h)/ ^JttY, necessary to observe 
the paramagnetic to ferromagnetic phase transition. This 
is plotted in Fig llOl versus the inverse exchange energy 
(1/Xoo). Here, 1/X^ is calculated from the electron 
density using Ea. (|ll|) . The corresponding r s values are 
also indicated on the top axis. The predicted dependence 
of D ss (Ep) on 1/Xoo, calculated using Ea. ljT^|> is plotted 
in Fie llUl for comparison. The agreement with theory is 
reasonably good, although there is some deviation at low 
electron density (high r s ). In these samples, spin split- 
ting occurs too early, before the critical density of states 
is reached. This may be a correction due to the non zero 
Zeeman energy. Clearly, in the limit of infinitely narrow 
Landau levels, our model breaks down, and the Zeeman 
energy is sufficient for the system to form a polarized 
ground state in the highest Landau level, even in the ab- 
sence of exchange interactions. The density of states for 
a 2D system in zero magnetic field is indicted by the dot- 
ted line in Fig^| The transition to a ferromagnetic state 
at zero magnetic field is predicted to occur for r s ~ 26, 
corresponding to the divergence of (x/xo)j m the calcu- 
lations of Attaccalite et al. 16 However, the corresponding 
electronic density of n s ~ 5 x 10 8 cm~ 2 is well below the 
density for the metal-insulator transition in GaAs^i At 
such densities, the system can no longer be described as 
a paramagnetic Fermi liquid, and hence no Stoner tran- 
sition is expected at B = in GaAs, in agreement with 
experiment. 

VI. CONCLUSION 

In conclusion, we have performed a quantitative anal- 
ysis of spin splitting in a series of GaAs heterojunctions 
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